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1 Introduction 

A significant part of today's researcli in 
QCD aims to provide tools to help bet- 
ter constrain the standard model and find 
what may lie beyond it. For example one 
wishes to determine, as accurately as possi- 
ble, the fundamentals of the QCD and elec- 
troweak theories, such as as, quark masses, 
and the elements of the CKM matrix. One 
also needs precise information about 'pseudo- 
fundamentals,' quantities such as parton dis- 
tribution functions (PDFs) that could be 
predicted if we knew how to solve non- 
perturbative QCD, but which currently must 
be deduced from experimental data. Finally, 
one puts this information together to predict 
the QCD aspects of both backgrounds and 
signals at high energy colliders, particularly 
at the Tevatron and LHC, to help maximise 
the chance of discovering and understanding 
any new physics. 

Other facets of QCD research seek to ex- 
tend the boundaries of our knowledge of QCD 
itself. The underlying Yang-Mills field theory 
is rich in its own right, and unexpected new 
perturbative structures have emerged in the 
past two years from considerations of string 
theory. In the high-energy limit of QCD it is 
believed that a new state appears, the widely 
studied colour glass condensate, which still 
remains to be well understood. And perhaps 
the most challenging problem of QCD is that 
of how to relate the partonic and hadronic 
degrees of freedom. 



Given the practical importance of QCD 
for the upcoming LHC programme, this talk 
will concentrate on results (mostly since the 
2003 Lepton-Photon symposium) that bring 
us closer to the well-defined goals mentioned 
in the first paragraph. Some of the more ex- 
plorative aspects will also be encountered as 
we go along, and one should remember that 
there is constant cross-talk between the two. 
For example: improved understanding of field 
theory helps us make better predictions for 
multi-jet events, which are important back- 
grounds to new physics; and by comparing 
data to accurate perturbative predictions one 
can attempt to isolate and better understand 
the parton-hadron interface. 

The first part of this writeup will be 
devoted to results at fixed order. At tree 
level we will examine new calculational meth- 
ods that are much more efficient than Feyn- 
man graphs; we will then consider NLO and 
NNLO calculations and look at the issues 
that arise in going from the Feynman graphs 
to useful predictions. 

One of the main uses of fixed-order order 
predictions is for understanding rare events, 
those with extra jets. In the second part of 
the writeup we shall instead turn to resum- 
mations, which help us understand the prop- 
erties of typical events. 

Throughout, the emphasis will be on the 
conceptual advances rather than the detailed 
phenomenology. Due to lack of space, some 
active current topics will not be covered, in 



particular exclusive QCD. Others are dis- 
cussed elsewhere in these proceedings.^ 

2 Fixed-order calculations 

2.1 Tree-level amplitudes and twistors 

Many searches for new physics involve signa- 
tures with a large number of final-state jets. 
Even for as basic a process as tt production, 
the most common decay channel (branching 
ratio of 46%), tt bbW+W^ -> bbqqqq in- 
volves 6 final-state jets, to which there are 
large QCD multi-jet backgrounds.^ And at 
the LHC, with 10fb~^ (1 year) of data, one 
expects of the order of 2000 events with 8 or 
more jets^ (pt(jet) > 60 GeV, % > 30deg, 

\y^\ < 3). 

For configurations with such large num- 
bers of jets, even tree-level calculations be- 
come a challenge — for instance gg — > 
8g involves 10525900 Feynman diagrams 
(see ref.^). In the 1980's, techniques 
were developed to reduce the complexity 
of such calculations.^ Among them colour 
decomposition,^ where one separates the 
colour and Lorentz structure of the ampli- 
tude, 

^'-^(l,2,...,n)=g"-2 ^ 

perms 

Tr(rir2...T„)A'-^(l,2,...,n); (1) 
\,^^^^^^_^ ^^^^^^^^ \,^^^^^^^_^^ ^^^^^^^^^ 

colour struct. colour ordered amp. 

the use^ of spinor products (ij) = = 
U-{ki)u+{kj) and [ij] = as the 

key building blocks for writing amplitudes; 
and the discovery^ and subsequent proof^ of 
simple expressions for the subset of ampli- 
tudes involving the maximal number of same- 
helicity spinors (maximum helicity violating 
or MHV amplitudes), i.e. n — 2 positive he- 
licity spinors for an n-gluon amplitude: 

Ref.^ also provided a computationally effi- 
cient recursion relation for calculating am- 



plitudes with arbitrary numbers of legs, and 
with these and further techniques,^" numer- 
ous programs (e.g. MadEvent," ALPGEN,^^ 
HELAC/PHEGAS,i3 

CompHEP,i4 GRACE, 15 Amegic^^) are able 
to provide results for processes with up to 10 
legs. 

The past two years have seen substantial 
unexpected progress in the understanding of 
multi-leg tree-level processes. It was initiated 
by the observation (first made by Nair^^) that 
helicity amplitudes have a particularly sim- 
ple form in 'twistor' space, a space where a 
Fourier transform has been carried out with 
respect to just positive helicity spinors. In 
twistor space a duality appears^^ between 
the weakly-coupled regimes of a topological 
string theory and TV = 4 SUSY Yang-Mills. 
This has led to the postulation, by Cac- 
hazo, Svrcek and Witten (CSW),^^ of rules 
for deriving non-MHV TV = 4 SUSY ampli- 
tudes from MHV ones, illustrated in fig. 1. 
For purely gluonic amplitudes the results are 
identical to plain QCD,^° because tree level 
SUSY amplitudes whose external legs are glu- 
ons have only gluonic propagators. 
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Figure 1. Graphical illustration of the CSW rules: 
by joining together two MHV amplitudes with an off- 
shell scalar propagator one obtains an amplitude with 
an extra negative helicity (NMHV). 

Recently a perhaps even more powerful 
set of recursion relations was proposed by 
Britto, Cachazo and Feng (BCF),^^ which 
allows one to build a general n-leg diagram 
by joining together pairs of on-shell sub- 
diagrams. This is made possible by contin- 
uing a pair of reference momenta into the 
complex plane. The proof^^ of these relations 
(and subsequently also of the CSW rules^'^) is 
remarkably simple, based just on the planar 



nature of colour-ordered tree diagrams, their 
analyticity structure, and the asymptotic be- 
haviour of known MHV ampHtudes. 

The discovery of the CSW and BCF 
recursion relations has spurred intense ac- 
tivity, about 150 articles citing the origi- 
nal papers^*'^*' having appeared in the 18 
months following their publication. Ques- 
tions addressed include the derivation of sim- 
ple expressions for specific amplitudes, the 
search for computationally efficient recursive 
formulations,'^''' extensions to processes with 
fermions and gluinos,^^ Higgs bosons,^'' elec- 
troweak bosons,^* gravity,^^ and the study of 
multi-gluon coUinear limits. This list is nec- 
essarily incomplete and further references can 
be found in a recent review^^ as well as below, 
when we discuss applications to loop ampli- 
tudes. 

2.2 One-loop amplitudes 

For quantitatively reliable predictions of 
a given process it is necessary for it to 
have been calculated to next-to-leading or- 
der (NLO). A wide variety of NLO calcu- 
lations exists, usually in the form of pub- 
licly available programs, that allow one to 
make predictions for arbitrary observables 
within a given process. The broadest of these 
programs are the MCFM,^^ NLOJET^^ and 
PHOX^s families. 

For a 2 ^ n process the NLO calcula- 
tion involves the 2 n + I tree-level dia- 
gram, the 2 ^ n 1-loop diagram, and some 
method for combining the tree-level matrix 
element with the loop contributions, so as to 
cancel the infrared and collinear divergences 
present in both with opposite signs. We 
have seen above that tree-level calculations 
are well understood, and dipole subtraction"^^ 
provides a general prescription for combining 
them with the corrresponding 1-loop contri- 
butions. The bottleneck in such calculations 
remains the determination of the 1-loop con- 
tribution. Currently, 2^3 processes are 



feasible, though still difficult, while as yet no 
full 2^4 1-loop QCD calculation has been 
completed. 

In view of the difficulty of these loop cal- 
culations, a welcome development has been 
the compilation, by theorists and experi- 
menters at the Los Houches 2005 workshop, 
Physics at TeV Colliders, of a realistic pri- 
oritised wish-list of processes. Among the 
most interesting remaining 2^3 processes 
one has pp WW+]ct, pp VVV {V = W 
or Z) and pp — > iJ -|- 2jets. The latter can 
be considered a 'background' to Higgs pro- 
duction via vector boson fusion, insofar as 
the isolation of the vector-boson fusion chan- 
nel for Higgs production would allow rela- 
tively accurate measurements of the Higgs 
couplings. A number of 2 — »■ 4 processes 
are listed as backgrounds to ttH production 
{pp ttqq,ttbb), WW H ^ WW or 
to general new physics {pp V + 3 jets) or 
specifically SUSY {pp VVV + jet). 

Two broad classes of techniques have 
been used in the past for 1-loop calcula- 
tions: those based directly on the evaluation 
of the Fc^ynman diagrams (sometimes for the 
1-loop-tree interference^^) and those based on 
unitarity techniques to sew together tree di- 
agrams (sec the rcvicw^°). Both approaches 
are still being actively pursued. 

Today's direct evaluations of 1-loop con- 
tributions have, as a starting point, the auto- 
mated generation of the full set of Feynman 
diagrams, using tools such as QGRAF^^ and 
FeynArts.*^ The results can be expressed in 
terms of sums of products of group-theoretic 
(e.g. colour) factors and tensor one-loop inte- 
grals, such as 

(3) 

Reduction procedures exist (e.g.'*'^''''') that 
can be applied recursively so as to express 

the In-fii..fj.i in terms of known scalar inte- 
grals. Such techniques form the basis of re- 
cent proposals for automating the evaluation 



of the integrals, where the recursion rela- 
tions are solved by a combination of ana- 
lytical and numerical methods,*^ or purely 
numerically, in some cases with spe- 
cial care as regards divergences that appear 
in the coefficients of individual terms of the 
recursion relation but vanish in the sum.'"''^^ 
Results from these approaches include a new 
compact form for gg 577 at 1 loop^*' 
and the 1-loop contribution to the 'priority' 
pp ^ H + 2 jets process. ^-'^ Related auto- 
mated approaches have also been developed 
in the context of electroweak calculations,^^ 
where recently a first 2^4 1-loop result 
{e~^e~ 4 fermions) was obtained. ^"^ 

Another approach^^ proposes subtrac- 
tion terms for arbitrary 1-loop graphs, such 
that the remaining part of the loop integral 
can carried out numerically in 4 dimensions. 
The subtraction terms themselves can be in- 
tegrated analytically for the sum over graphs 
and reproduce all infrared, coUinear and ul- 
traviolet divergences. 

The above methods are all subject to the 
problem of the rapidly increasing number of 
graphs for multi-leg processes. On the other 
hand, procedures based on 'sewing' together 
tree graphs to obtain loop graphs^° (cut con- 
structibility approach) can potentially bene- 
fit from the simplifications that emerge from 
twistor developments for tree graphs. This 
works best for TV = 4 SUSY QCD, where 
cancellations between scalars, fermions and 
vector particles makes the 'sewing' procedure 
simplest and for example all gluonic (and 
some scalar and gluino) NMHV 1-loop he- 
licity amplitudes are now known;^^ also, con- 
jectures for TV = 4 SUSY n-leg MHV pla- 
nar graphs, at any number of loops, based 
on 4-gluon two and three-loop calculations,"''^ 
have now been explicitly verified for 5 and 
6 gluon two-loop amplitudes. For Af = 
1 SUSY QCD, known results for all MHV 
amplitudes'^ have been reproduced'^ and 
new results exist for some all-n NMHV 
graphs^° as well as full results for 6 gluons.^^ 



For plain QCD, progress has been slower, 
though all finite 1-loop graphs (-|- -|- -|- -|- . . ., 
— !- + + •••) were recently presented^^ and 
understanding has also been achieved for di- 
vergent graphs, including the full result for 
all 1-loop (- + . . .) MHV graphs.6^ The 
prospects for the twistor-inspired approach 
are promising and one can hope that it will 
soon become practically competitive with the 
direct evaluation methods. 

2.3 NNLO jet calculations 

NNLO predictions are of interest for many 
reasons: in those processes where the pertur- 

bative series has good convergence they can 
help bring perturbative QCD predictions to 
the percent accuracy level. In cases where 
there are signs of poor convergence at NLO 
they will hopefully improve the robustness 
of predictions and, in all cases, give indica- 
tions of the reliability of the series expan- 
sion. Finally they may provide insight in 
the discussion of the relative importance of 
hadronisation and higher-order perturbative 
corrections. ^''^^'^'^ 

So far, NNLO predictions are available 
mostly only for processes with 3 external 
legs, such as the total cross section for Z 
hadrons or the inclusive pp W, ^68,69,70 ^j. 
jjjggg69,7i,72 pj.Qgg sections. A full list is given 
in table 1 of Stirling's ICHEP writeup.^^ 

Most current effort is being directed to 
the e+e^ 3jets process, where NLO cor- 
rections are often large, and where one is 
free of complications from incoming coloured 
particles. The ingredients that are needed 
are the squared 5-parton tree level (M5) and 
3-parton 1-loop (M30) amplitudes, and the 
interference between 4-parton tree and 1- 
loop (M4), and between 3-parton tree and 2- 
loop amplitudes (M^i,), all of which are now 
known (for references see introduction of^^). 
A full NNLO prediction adds the integrals 
over phase-space of these contributions, mul- 
tiplied by some jet-observable function J that 



depends on the momenta Pi 
Jnnlo = j dP^^Mr,J{pi„^) + 

J d^$4M4 J(pi..4)+ J d^$3(M3a+6) J(P1..3) • 

Each of the terms is infrared and coUincar 
(IRC) divergent, because of the phase space 
integration (M4_5) and/or the loop integral 
(-^3,4)- The current bottleneck for such cal- 
culations is in cancelling these divergences for 
an arbitrary (IRC safe) J. 

A standard approach at NLO is to intro- 
duce subtraction terms (e.g. in the dipole^^ 
formalism) , schematically, 

+ J A(M4J(pi..4)+54J(pi..4)) 

such that each integral is separately finite 
and that, alone, the ^5 and 5*4 terms inte- 
grate to equal and opposite divergent con- 
tributions (they both multiply the same 4- 
parton jet function J{pi..4) and the P1..4 are 
a specifically designed function of the P1..5). 
At NNLO a similar method can be envisaged, 
and considerable work has gone towards de- 
veloping a general formalism. 

An alternative 
approach, sector decomposition, in- 
troduces special distributions /_j (involving 
plus-functions, like those in splitting func- 
tions), which isolate the e~* divergent piece 
of a given integral (e = (D — 4)/2), 

J d%M5J(pi..5) = ^J d%5/-4M5J(/)l..5) 
+ d%5/-3M5J(pi..5) + ... , 

where the integration is performed in trans- 
formed variables that simplify the separation 

of divergences. In such an approach one thus 
obtains separate results for each power of e in 
both real and virtual terms, making it easy 
to combine them. 

A useful testing ground for a num- 
ber of these approaches has been e+e^ 



2 jets. Fundamentally new results are 
the differential distributions for W, Z, H 
production^^ in the sector-decomposition ap- 
proach and the (a^ (7^/277)^ contribution to 
(1 - Thrust) in e+e" 3 jets, -20.4 ± 4, in 
the 'antenna' subtraction approach.^^ 

2.4 NNLO splitting functions 

The landmark calculation of 2004 was prob- 
ably that of the NNLO splitting functions by 
Moch, Vermaseren and Vogt^'^ (MVV). These 
are important for accurate DGLAP evolution 
of the parton distributions, as extracted from 
fixed target, HERA and Tevatron data, up to 
LHC scales. 

The results for the splitting functions 
take about 15 pages to express, though the 
authors have also provided compact approx- 
imations for practical use. These are gradu- 
ally being adopted in NNLO fits.*^'^^*'^ Mostly 
the NNLO splitting functions are quite sim- 
ilar to the estimates obtained a few years 
ago^^ based on a subset of the moments and 
known asymptotic limits. In particular it re- 
mains true that the NNLO corrections are in 
general small, both compared to NLO and in 
absolute terms. The worst region is that of 
small X for the singlet distributions, shown in 
fig. 2, where the NLO corrections were large 
and there is a significant NNLO modification 
as well. Studies of small-x resummation sug- 
gest that further high-order effects should be 
modest, ^'''^^ though so far only the gluon sec- 
tor has been studied in detail. 

Another potentially dangerous region is 
that of large x: the splitting functions con- 
verge well, but the coefficient functions have 
[a"ln^"'~^(l — x)]/{l — x) enhancements. 
There are suggestions that the all-order in- 
chision of these enhancements via a thresh- 
old resummation may help improve the accu- 
racy of PDF determinations.^^ Fresh results 
from the MVV group, ^"^ for the third order 
electromagnetic coefficient functions, quark 
and gluon form factors, and new threshold 
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Figure 2. Impact*^ of NNLO DGLAP corrections on 
the derivative of a toy singlet quark distribution qs = 

resummation coefficients should provide the 
necessary elements for yet higher accuracy at 
large x. Together with the large-x part of 
the NNLO splitting functions these have been 
used as inputs to the N'^LO soft-collinear en- 
hanced terms for the Drell-Yan and Higgs 
cross sections. ^^'^^ 

Finally, as splitting-function calcula- 
tions approach accuracies the 1% level, one 
should consider also the relevance of QED 
corrections. ^^'^^ 

2.5 Other accuracy-related issues 

In view of the efforts being devoted to im- 
proving the accuracy of theoretical predic- 
tions, it is disheartening to discover that 
there are still situations where accuracy 
is needlessly squandered through incorrect 
data-theory comparisons. This is the case 
for the inclusive jet cross sections in the cone 
algorithm at the Tevatron, where a parame- 
ter Rsep ~ 1.3 is introduced in the cone al- 
gorithm used for the NLO calculation (see 
e.g.^^), but not in the cone algorithm ap- 
plied to the data. Rsep is the multiple of the 
cone radius beyond which a pair of partons is 
not recombined. It dates to early theoretical 



work^^ on NLO corrections to the cone al- 
gorithm in hadron collisions, when there was 
no public information on the exact jet algo- 
rithm used by the experiments: Rsep was in- 
troduced to parametrise that ignorance. 

The use of Rsep in just the theory in- 
troduces a spurious NLO correction (at the 
5 — 10% level^^), meaning that the data- 
theory comparisons are only good to LO. 
The size of the discrepancy is comparable 
to the NLO theory uncertainty. As this is 
smaller than experimental errors, for now the 
practical impact is limited. However, as ac- 
curacies improve it is essential that theory- 
experiment comparisons be done consistently, 
be it with properly used cone algorithms'^ or 
with the (more straightforward and power- 
ful) kt algorithm, ^"^^ which is finally starting 
to be investigated.^"^ 

An accuracy issue that is easily over- 
looked when discussing QCD developments 
is the non-negligible impact of electroweak ef- 
fects at large scales. The subject has mostly 
been investigated for leptonic initial and final 
states at a linear collider (ILC). The dom- 
inant contributions go as asw^^^ Pt/Mw- 
Since the LHC can reach transverse momenta 
(Pf) an order of magnitude larger than the 
ILC the electroweak effects are very consid- 
erably enhanced at the LHC, being up to 
30 - 40%}°^'^^^ Among the issues stiU to be 
understood in such calculations (related to 
the cancellation of real and virtual correc- 
tions) is the question of whether experiments 
will include events with W and Z^s as part of 
their normal QCD event sample, or whether 
instead such events will be treated separately. 

3 All-order calculations 

All-order calculations in QCD are based on 
the resummation of logarithmically dominant 
contributions at each order. Such calcula- 
tions are necessary if one is to investigate the 
properties of typical events, for which each 
extra power of is accompanied by large 



soft and coUinear logarithms. 

The two main ways of obtaining all- 
order resummed predictions are with exclu- 
sive Monte Carlo event generators, and with 
analytical resummations. The former provide 
moderately accurate (leading log (LL) and 
some parts of NLL) predictions very flexibly 
for a wide variety of processes, with hadroni- 
sation models included. The latter are able 
to provide the highest accuracy (at parton 
level), but usually need to be carried out 
by hand (painfully) for each new observable 
and/or process. All-order calculations are 
also used in small- a; and saturation physics, 
which will be discussed briefly at the end of 
this section. 

3.1 Monte Carlo event generators (MC) 

Various issues arc present in current work on 
event generators — the switch from Fortran 
to C-|— 1-; improvements in the showering al- 
gorithms and the modelling of the underlying 
event; and the inclusion of information from 
fixed-order calculations. 

The motivation for moving to C-|— I- is the 
need for a more modern and structured pro- 
gramming language than Fortran 77. C-|— |- 
is then a natural choice, in view both of its 
fl(^xibility and its widespread use in the ex- 
perimental community. 

Originally it was 

intended that Herwig^"'', Pythia^°^ and 
Ariadne/LDC^"*''^°'' should all make use of 
a general C-| — h event generator framework 
known as ThePEG.i°s Herwig-|-F,i°9 based 
on ThePEG, was recently released for e"'"e~ 
and work is in progress for a hadron-hadron 
version. Pythia 7 was supposed to have been 
the Pythia successor based on the ThePEG, 
however instead a standalone C++ gener- 
ator Pythia 8 is now being developed, 
perhaps to be interfaced to ThePEG later 
on. Another independent C++ event gen- 
erator has also recently become available, 
SHERPA,^^^ whose showering and hadroni- 



sation algorithms arc largely based on those 
of Pythia, and which is already functional for 
hadron collisions. 

In both the Pythia and Herwig 'camps' 
there have been developments on new show- 
ering algorithms. Herwig++ incorporates an 
improved angular-ordered shower^^^ in which 
the 'unpopulated' phase space regions have 
been shrunk. Pythia 6.3 has a new par- 
ton shower^^^ based on transverse momen- 
tum ordering (i.e. somewhat like Ariadne) 
which provides an improved description of 
e+e~ event-shape data and facilitates the 
modelling of multiple interactions in hadron 
collisions. Separately, investigations of al- 
ternatives to standard leading-log backward 
evolution algorithms for initial-state showers 
are also being pursued. ^^^'^^^ 

3.2 Matching MC & fixed order 

Event generators reproduce the emission pat- 
terns for soft and coUinear gluons and also 
incorporate good models of the transition to 
hadrons. They are less able to deal with 
multiple hard emissions, which, as discussed 
above, are important in many new particle 
searches. There is therefore a need to com- 
bine event generators with fixed-order calcu- 
lations. 

The main approach for this is the 
CKKW^^^ proposal. Events are generated 
based on the n-parton tree matrix elements 
(for various n), keeping an event only if its 
n partons are sufficiently well separated to 
be considered as individual jets (according to 
some threshold jet-distance measure, based 
e.g. on relative kt). Each event is then as- 
signed a 'best-guess' branching history, us- 
ing which it can be reweighted with appro- 
priate Sudakov form factors (to provide vir- 
tual corrections) and running couplings. The 
normal parton showering is then added on to 
each event at scales below the threshold jet- 
distance measure. Over the past two years 
this method has become widely adopted and 



is available in all major generators. 

As well as seeking to describe more jets 
it is also important to increase the accuracy 
of event generators for limited number of jets, 
by including NLO corrections. Here too there 
is one method that has so far dominated 
practical uses, known as MC@NLO.^^^ Very 
roughly, it takes the standard MC and mod- 
ifies it according to 

MC^MC(1 + NL0-NL0mc), (4) 

where NLOmc represents the effective NLO 
corrections present by default in the MC. 
Since the MC usually has the correct soft 
and coUinear divergences, the combination 
NLO - NLOmc should be finite. Eq. (4) 
is therefore a well-behaved way of introduc- 
ing exactly the correction needed to guar- 
antee NLO correctness. A prediction from 
MC@NLO for 6-production,"^ figure 3, is 
compared to data and a purely analytical 
approach, and one notes the good agree- 
ment between all three. 







- Points: CDF 




Curves: KONl.I, ^SSk 




; (j(pj(J/lI')>1.25 GeV) BR; ^^i^S^ 




: 19.9-3:2 nb (CDF) ^^^J^S^l^ 




18.3^53 nb (FONLL) 




: Solid iiiBtogram: MCigiNLO. 17.2 nb. 




- Dashed histogram: MCeNLO, 16.4 nb 





5 10 15 20 

p.,(J/li) (GeV) 

Figure 3. Spectrum for b-production {— + J/^) com- 
pared to MCQNLO and an analytical prediction. -"^^^ 

A bottleneck in widespread implementa- 
tion of the MC@NLO approach is the need 
to know NLOmc J which is different for each 
generator, and even process. Furthermore 
MC@NLO so far guarantees NLO correct- 
ness only for a fixed number of jets — e.g. 
it can provide NLO corrections to W pro- 
duction, but then W + 1 jet is only provided 



to LO. An approach to alleviate both these 
problems proposes^^^ to combine, for exam- 
ple, W, M^-l- 1 jet, VK-t-2jets, etc. with a pro- 
cedure akin to CKKW. It seeks to alleviate 
the problem of needing to calculate NLOmc 
as follows: when considering W + mjets the 
(m-t- 1)**^ emission (that needed for NLO ac- 
curacy) is generated not by the main MC, 
but by a separate mini well-controlled gen- 
erator, designed specifically for that purpose 
and whose NLO expansion is easily calcu- 
lated (as is then the analogue of eq.(4)). The 
only implementation of this so far (actually 
of an earlier, related formalism^^^) has been 
for e+e-.i23 

3.3 Analytical resummations 

It is in the context of analytical resumma- 
tions that one can envisage the highest re- 
summation accuracies, as well as the sim- 
plest matching to fixed order calculations. 
Rather than directly calculating the distribu- 
tion d(j(y)/dV of an observable V , one often 
considers some integral transform F{v) of the 
distribution so as to reduce F{v) to the form 

lnF(:.) = ^(<_L2+^+<L^+...), (5) 

" LL NLL 

with L = \nv. 

Much of the information for certain 
N'^LL threshold resummations was recently 
provided by the MVV group. The high- 
est accuracy for full phenomenological dis- 
tributions is for the recently calculated 
Higgs transverse momentum spectrum^^'* at 
the LHC and the related^^^ energy-energy- 
correlation (EEC) in e"'"e~,^^^ both of which 
have also been matched to NLO fixed order. 

Boson transverse momentum spectra and 
the EEC are among the simplest observables 
to resum. For more general observables and 
processes, such as event shapes or multi-jet 
events, the highest accuracy obtained so far 
has been NLL, and the calculations are both 
tedious and error-prone. This has prompted 



work on understanding general features of re- 
summation. One line of research^^^ exam- 
ines the problem of so-called 'factorisable' ob- 
servables, for an arbitrary process. This ex- 
tends the understanding of large-angle soft 
colour evolution logarithms, whose resumma- 
tion was originally pioneered for 4-jet pro- 
cesses by the Stony Brook group, ^^^'^^^ in 
which unexplained hidden symmetries were 
recently discovered, "'^'^"'^^^ notably between 
kinematic and colour variables. 

Separately, the question of how to treat 
general observables has led to a procedure for 
automating resummations for a large class 
of event-shape- like observables. ^'^^ It avoids 
the need to find an integral transform that 
factorises the observable and introduces a 
new concept, recursive infrared and coUinear 
safety, which is a sufficient condition for the 
exponentiated form eq. (5) to hold. Its main 
application so far has been to hadron-coUider 
dijet event shapes, ^^'^ (see also^'^^) which pro- 
vide opportunities for experimental investiga- 
tion of soft-colour evolution and of hadronisa- 
tion and the underlying event at the Tevatron 
and LHC. 

The above resummations apply to 
'global' observables, those sensitive to radi- 
ation everywhere in an event. For non-global 
observables, such as gap probabilities^^^ or 
properties of individual jets,^^^'^"^^ a new 
class of enhanced term appears, non-global 
logarithms a"L" (NGL), fig. 4. Their resum- 
mation has so far only been possible in the 
large- A^c limit, ^^^'^'^^ though the observation 
of structure related to BFKL evolution, ^'^^ 
has inspired proposals for going to finite 

Non-global and soft colour resumma- 
tions (in non-inclusive form^^^) come to- 
gether when calculating the probability of 
a gap between a pair of jets at the Teva- 
tron/LHC, relevant as a background to the 
W fusion process for Higgs production. ^''^ 
Recently it has been pointed out that contri- 
butions that are subleading in play an 




Figure 4. Diagram giving NGL: to calculate the prob- 
ability of there being no emission into the gap, one 
should resum a large-angle energy-ordered cascade of 
emissions rather than just direct emission from the 
original hard partons. 

important role for large gaps,^"*^ and also that 
there are considerable subtleties when using 
a kt jet algorithm to to define the gap.-^**^ 

3.4 Small X and saturation 

The rise of the gluon at small a;, as pre- 
dicted by BEKL,^"*^ leads eventually to such 
high gluon densities that a 'saturation' phe- 
nomenon should at some point set in. Usu- 
ally one discusses this in terms of a satura- 
tion scale, Ql{x), below (above) which the 
gluon distribution is (un)saturated. From 
HERA data, it is believed that at a; = xq ~ 
10""* - 10~^ Qs is of order 1 GeV and that it 
grows as Ql{x) = {x/xo)~'^ GeV^ with A ~ 
0.3.^**^ Qs niay be of relevance to the LHC 
because the typical transverse scale Em of 
minimum bias minijets should satisfy the re- 
lation Q1{s/E'^) ~ whose solution gives 
Em ~ (sa;o)5TA GeV^^ ~ 2.7 - 3.6 GeV, or 
if one doesn't trust the normalisation, a fac- 
tor 1.7 relative to the Tevatron. The exact 
phenomenology is however delicate. 

The theoretical study of the saturation 
scale has seen intense activity these past 18 
months, spurred by two observations. Firstly 
it was pointed out^''^ that the Balitsky- 
Kovchcgov (BK) equation, often used to de- 
scribe the onset of saturation and the evolu- 
tion dlnQg/dlnx, is in the same universal- 
ity class as the Fisher Kolmogorov Petrovsky 
Piscounov (FKPP) reaction-diffusion equa- 
tion, much studied in statistical physics""^*^ 
and whose travelling wave solutions relate to 



the evolution of the saturation scale. Sec- 
ondly, large corrections were discovered^^° 
when going beyond the BK mean-field ap- 
proximation: to LO, Xbk oc as, while the 
non mean-field corrections go as Us/ln^al. 
Such corrections turn out to be a famil- 
iar phenomenon in stochastic versions of the 
FKPP equation, with in QCD playing 
the role of the minimum particle density in 
reaction-diffusion systems with a finite num- 
ber of particles. ^■''^ The stochastic corrections 
also lead to a large event-by-event dispersion 
in the saturation scale. 

These stochastic studies are mostly based 
on educated guesses as to the form of the 
small- a; evolution beyond the mean-field ap- 
proximation. There has also been extensive 
work on finding the full equation that re- 
places BK beyond the mean-field approxima- 
tion, a number of new formulations having 
been proposed.-^^^'^^^'-^^^'-^^^ It will be inter- 
esting to examine how their solutions com- 
pare to the statistical physics related ap- 
proaches. 

4 Concluding remarks 

Of the topics covered here, the one that has 
been the liveliest in the past year is that 
of 'twistors'." Its full impact cannot yet be 
gauged, but the dynamic interaction between 
the QCD and string-theory communities on 
this subject will hopefully bring further im- 
portant advances. 

More generally one can ask if QCD is on 
track for the LHC. Progress over the past 
years, both calculationally and phenomeno- 
logically (e.g. PDF fitting) has been steady. 
Remaining difficulties, for example in high 
(NNLO) and moderate (NLO) accuracy cal- 
culations are substantial, however the consid- 
erable number of novel ideas currently being 
discussed encourages one to believe that sig- 
nificant further advances will have been made 
by the time LHC turns on. 

"Second place goes to saturation. 
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DISCUSSION 

Bennie Ward (Baylor University): 

To what extent are the deduced MHV 
rules using twistors now proven? 

Gavin Salam: an outline of a field-theoretic 
proof for the MHV (CSW) rules was 
given in the same article^^ as the proof of 
the BCF rules. Recently (after the Lep- 
ton Photon Symposium) a more detailed 
version of the proof has appeared. 



